3. The Motion of Rigid Bodies 




Figure 22: Wolfgang Pauli and Niels Bohr stare in wonder at a spinning top. 



Having now mastered the technique of Lagrangians, this section will be one big 
application of the methods. The systems we will consider are the spinning motions of 
extended objects. As we shall see, these can often be counterintuitive. Certainly Pauli 
and Bohr found themselves amazed! 



We shall consider extended objects that don't have any internal 
degrees of freedom. These are called "rigid bodies" , defined to be 
a collection of N points constrained so that the distance between 
the points is fixed, i.e. 



r,-| = constant 



(3.1) 




Figure 23: 



for all i,j — l,...,N. A simple example is a dumbell (two masses 
connected by a light rod), or the pyramid drawn in the figure. In both cases, the 
distances between the masses is fixed. 
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Often we will work with continuous, rather than discrete, bodies simply by replacing 
J2i m i ~* I dr p(r) where p(r) is the density of the object. A rigid body has six degrees 
of freedom 

3 Translation + 3 Rotation 

The most general motion of a free rigid body is a translation plus a rotation about 
some point P. In this section we shall develop the techniques required to describe this 
motion. 

3.1 Kinematics 

Consider a body fixed at a point P. The most general allowed motion is a rotation 
about P. To describe this, we specify positions in a fixed space frame {e a } by embedding 
a moving body frame {e a } in the body so that {e a } moves with the body. 




e 3 (t 2 ) 



Figure 24: The fixed space frame and the moving body frame. 

Both axes are orthogonal, so we have 

e a • e 6 = 5 ab , e a (t) • e b (t) = 5 ab (3.2) 
We will soon see that there is a natural choice of the basis {e a } in the body. 

Claim: For all t, there exists a unique orthogonal matrix R(t) with components R a b{t) 
such that e a (t) = R ab {t)e b 

Proof: e a ■ e b = 5 ab =>- R ac R bd e c ■ e d = 5 ab => R ac Rbc = S ab or, in other words, 
(R T R) ab = 5 ab which is the statement that R is orthogonal. The uniqueness of R 
follows by construction: R ab = e a ■ e b . □. 
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So as the rigid body rotates it is described by a time dependent orthogonal 3x3 
matrix R(t). This matrix also has the property that its determinant is 1. (The other 
possibility is that its determinant is —1 which corresponds to a rotation and a reflection 
e a — > — e ). Conversely, every one-parameter family R(t) describes a possible motion 
of the body. We have 

C = Configuration Space = Space of 3 x 3 Special Orthogonal Matrices = 5*0(3) 

A 3 x 3 matrix has 9 components but the condition of orthogonality R T R = 1 imposes 
6 relations, so the configuration space C is 3 dimensional and we need 3 generalised 
coordinates to parameterise C. We shall describe a useful choice of coordinates, known 
as Euler angles, in section 3.5. 

3.1.1 Angular Velocity 

Any point r in the body can be expanded in either the space frame or the body frame: 

r (t) = ?a(t) e a in the space frame 

= r a e a (t) in the body frame (3.3) 

where fb(t) = r a R ab (t). Taking the time derivative, we have 

dr df a 

— = —— e a m the space frame 
at dt 

= r a — in the body frame 
dt J 

dRab _ (r> ,x 

= r a — e b (3.4) 

Alternatively, we can ask how the body frame basis itself changes with time, 

de a dR ab _ f dR ab A , N 

e 6 = ( — — R e c = Lo ac e c (3.5) 

/ be 



dt dt \ dt 



where, in the last equality, we have defined u ac = R a b(R 1 )fec = RabRcb using the fact 
that R T R = 1. 

Claim: u ac = —uj ca i.e. u is antisymmetric. 

Proof: R ab R cb = 5 ac =5- R ab R cb + R ab R cb = =3- uj ac + u ca = □ 
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Since uj ac is antisymmetric, we can use it to define an object with a single index 
(which we will also call a;) using the formula 



— 2 € abc^bc 



(3.6) 



so that UJ3 = ujyi and so on. We treat these u a as the components of a vector in the 
body frame, so u> = u a e a . Then finally we have our result for the change of the body 
frame basis with time 



dt 



-e abc uj b e c = uj x e 



(3.7) 



where, in the second equality, we have used the fact that our body frame axis has 
a "right-handed" orientation, meaning e a x e b = e abc e c . The vector u> is called the 
instantaneous angular velocity and its components ui a are measured with respect to the 
body frame. 

Since the above discussion was a little formal, let's draw a picture 
to uncover the physical meaning of u>. Consider a displacement of a 
given point r in the body by rotating an infinitesimal amount d(f> about 
an axis n. From the figure, we see that |dr| = |r| d(f) sin 9. Moreover, 
this displacement is perpendicular to r since the distance to P is fixed 
by the definition of a rigid body. So we have 

dr = d<j> x r (3.8) 

with d<f> = nd(f). "Dividing" this equation by dt, we have the result 

r = u) x r 




Figure 25: 



(3.9) 

where u> = d<p/dt is the instantaneous angular velocity. In general, both the axis of 
rotation n and the rate of rotation d(j>/dt will change over time. 

Aside: One could define a slightly different type of angular velocity by looking at 
how the space frame coordinates f a (t) change with time, rather than the body frame 
axes e a . Since we have r b (t) = r a R ab (t), performing the same steps as above, we have 



r b = r a R ba = r a (R R) ab 



(3.10) 



which tempts us to define a different type of angular velocity, sometimes referred to as 
"convective angular velocity" by Q ab = R~^R cb which has the and R in a different 
order. Throughout our discussion of rigid body motion, we will only deal with the 
original 00 = RR^ 1 . 



- 48 - 



3.1.2 Path Ordered Exponentials 

In the remainder of this chapter, we will devote much effort to determine the the angular 
velocity vector u>(t) of various objects as they spin and turn. But how do we go from 
this to the rotation R(t)7 As described above, we first turn the vector u) = w a e a into a 
3x3 antisymmetric matrix u a b = e a bcU c - Then, from this, we get the rotation matrix 
R by solving the differential equation 

u = ^R- 1 (3.11) 
dt y ' 

If oj and R were scalar functions of time, then we could simply integrate this equation 
to get the solution 

R(t) = exp ^ jf uj(t') dt'^j (3.12) 

which satisfies the initial condition R(0) = 1. But things are more complicated because 
both to and R are matrices. Let's first describe how we take the exponential of a matrix. 
This is defined by the Taylor expansion. For any matrix M, we have 

exp(M) = 1 + M + ±M 2 + ... (3.13) 

As our first guess for the solution to the matrix equation (3.11), we could try the scalar 
solution (3.12) and look at what goes wrong. If we take the time derivative of the 
various terms in the Taylor expansion of this putative solution, then problems first 
arise when we hit the \M 2 type term. The time derivative of this reads 

\ jt Of u{t>) dt ) = \ u ® Of dt ) + \ dt ) ^ (3,i4) 

We would like to show that R = uoR. The first term on the right-hand side looks 
good since it appears in the Taylor expansion of ujR. But the second term isn't right. 
The problem is that we cannot commute u(t) past cu(t') when t' ^ t. For this reason, 
equation (3.12) is not the solution to (3.11) when oo and R are matrices. But it does 
give us a hint about how we should proceed. Since the problem is in the ordering of 
the matrices, the correct solution to (3.11) takes a similar form as (3.12), but with a 
different ordering. It is the path ordered exponential, 

R(t) = P exp ^ jf uj(t') dt'^j (3.15) 
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where the P in front means that when we Taylor expand the exponential, all matrices 
are ordered so that later times appear on the left. In other words 

R(t) = l+ / u(t')dt' + / / uj{t")uj{t') dt'dt" + ... (3.16) 
Jo Jo Jt' 

The double integral is taken over the range < t' < t" < t. If we now differentiate 
this double integral with respect to t, we get just the one term iv(t) ^ f Q * u(t') dt'^j , 
instead of the two that appear in (3.14). It can be checked that the higher terms in the 
Taylor expansion also have the correct property if they are ordered so that matrices 
evaluated at later times appear to the left in the integrals. This type of path ordered 
integral comes up frequently in theories involving non-commuting matrices, including 
the standard model of particle physics. 

As an aside, the rotation matrix R is a member of the Lie group SO (3), the space of 
3x3 orthogonal matrices with unit determinant. The antisymmetric angular velocity 
matrix u, corresponding to an instantaneous, infinitesimal rotation, lives in the Lie 
algebra so(3). 

3.2 The Inertia Tensor 

Let's look at the kinetic energy for a rotating body. We can write 

i 

= ^m 1 (uxr 1 )-(wx n) 

i 

= iX>((w-aOfc-*)-fo-«) 2 ) ( 3 - 17 ) 

i 

Or, in other words, we can write the kinetic energy of a rotating body as 

T = \u a I ab uj b (3.18) 

where I a b, a,b — 1, 2, 3 are the components of the inertia tensor measured in the body 
frame, defined by 

Lb = ^2 m i (( r i • r i)$ab - (ri)o(ri)b) (3.19) 

i 

Note that I a b = h a s° the inertia tensor is symmetric. Moreover, the components 
are independent of time since they are measured with respect to the body frame. For 
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continuous bodies, we have the analogous expression 



J dh p(r) 



( y 2 + z 2 



-xy 



—xy x 2 + z 2 



-xz 
-yz 

2 i „,2 



\ 



(3.20) 



y — xz —yz x + y J 



Since I a b is a symmetric real matrix, we can diagonalise it. This means that there 
exists an orthogonal matrix O such that OIO T = I' where I' is diagonal. Equivalently, 
we can rotate the body frame axis {e a } to coincide with the eigenvectors of / (which 
are {Oe a }) so that, in this frame, the inertia tensor is diagonal. These preferred body 
axes, in which / is diagonal, are called the principal axes. In this basis, 



h 



\ 



\ 



(3.21) 



The eigenvalues I a are called the principal moments of inertia. The kinematical prop- 
erties of a rigid body are fully determined by its mass, principal axes, and moments of 
inertia. Often the principal axes are obvious by symmetry. 



Claim: The I a are real and positive. 
Proof: If c is an arbitrary vector, then 



I ab c a c b = ™i(r-c 2 - (Ti ■ c) 2 ) > (3.22) 

i 

with equality only if all the lie on a line. If c is the a th eigenvector of I then this 
result becomes I a bC a c b = I a \c\ 2 which tells us I a > 0. □ 

Example: The Rod 

Consider the inertia tensor of a uniform rod of length / and mass M about its centre. 
The density of the rod is p = M/l. By symmetry, we have / = diag(/i, ii, 0) where 

nl/2 l 

h = I px 2 dx= —Ml 2 (3.23) 
J -i/2 12 

Example: The Disc 
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Now consider a uniform disc of radius r and mass M. 
We take the z axis to be perpendicular to the disc and mea- 
sure / about its centre of mass. Again we know that / = 
diag(7i, J 2 , h)- The density of the disc is p = M/irr 2 , so we 
have 



h = / py 2 d 2 x , I 



-s 



2 j2 

px a x 




Figure 26: 



so ii = I 2 by symmetry, while 



h = J p(x 2 + y 2 )d 2 x 



Therefore 



I 3 = I l + I 2 = 2np I r' 3 dr' = \Mr 2 



(3.24) 



So the moments of inertia are I\ — I 2 — \Mr 2 and I 3 = \Mr 2 . 
3.2.1 Parallel Axis Theorem 

The inertia tensor depends on what point P in the body is held fixed. In general, if we 
know I about a point P it is messy to compute it about some other point P'. But it 
is very simple if P happens to coincide with the centre of mass of the object. 

Claim: If P' is displaced by c from the centre of mass, then 



{Qab = (I C .of.m)ab + M(c 2 5 ab - C a C b ) 



(3.25) 



Proof: 



(Ic)ab = mi i ( Ti ~ C ^ 5ah ~ ( Ti ~ ~~ ( 3 - 26 ) 

i 

= ^2m {riS ab - (r;) a (r;) 6 + [-2Ti ■ c5 ab + (ri) a c 6 + (rj) 6 c a ] + (c 2 5 ab - c a c 6 )} 



But the terms in square brackets that are linear in vanish if rj is measured from the 
centre of mass since £V m^i = 0. □ 

The term M(c 2 5 ab — c a c b ) is the inertia tensor we would find if the whole body was 
concentrated at the centre of mass. 
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Example: The Rod Again 

The inertia tensor of the rod about one of its ends is l\ = -^Ml 2 + M(//2) 2 = |MZ 2 . 
Example: The Disc Again 

Consider measuring the inertia tensor of the disc about a point displaced by c = (c, 0, 0) 
from the centre. We have 



/ l r 2 



M 



1„2 
4' 



1„2 
2' 



+ M 



( ±r 2 



= M 



\ 




l r 2 + c 2 



V 



\r 2 + c 2 



Figure 27: 



3.2.2 Angular Momentum 

The angular momentum L about a point P can also be described neatly in terms of 
the inertia tensor. We have 

i 

= ^m^rj x (u; x r,) 

i 

= ^rriiirfa - (u> ■ TjTi) 

i 

= Iu> (3.27) 

In the body frame, we can write L = L a e a to get 

L a = I ab u b (3.28) 

where uj = oo a e a . Note that in general, u> is not equal to L: the spin of the body and 
its angular momentum point in different directions. This fact will lead to many of the 
peculiar properties of spinning objects. 

3.3 Euler's Equations 

So far we have been discussing the rotation of a body fixed at a point P. In this section 
we will be interested in the rotation of a free body suspended in space - for example, a 
satellite or the planets. Thankfully, this problem is identical to that of an object fixed 
at a point. Let's show why this is the case and then go on to analyse the motion. 



- 53 - 



Figure 28: 

The most general motion of a body is an overall translation superposed with a ro- 
tation. We could take this rotation to be about any point in the body (or, indeed, a 
point outside the body). But it is useful to consider the rotation to be about the center 
of mass. We can write the position of a particle in the body as 

Ti(t)=R(t) + ATi(t) (3.29) 

where Arj is the position measured from the centre of mass. Then examining the 
kinetic energy (which, for a free body, is all there is) 

T = \Y j m l i 2 l 

i 

= Yl mi + R • x Ar i) + >x Ar,) 2 

i 

= \MYC + \u a I ab u> b (3.30) 

where we've used the fact that £V mjAi-j = 0. So we find that the dynamics separates 
into the motion of the centre of mass R, together with rotation about the centre of 
mass. This is the reason that the analysis of the last section is valid for a free object 

3.3.1 Euler's Equations 

From now on, we shall neglect the center of mass and concentrate on the rotation of 
the rigid body. Since the body is free, its angular momentum must be conserved. This 
gives us the vector equation 

f = (3,!) 
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Let's expand this in the body frame, we have 

Q = dL = dLa_ e | L de a 

dt dt a a dt 

= ^e a + L a u:xe a (3.32) 

This simplifies if we choose the body axes {e a } to coincide with the the principal axes. 
Using L a = I a bUJb, we can then write L x = 1\UJ\ and so on. The equations of motion 
(3.32) are now three non-linear coupled first order differential equations, 

Iiiui + Ld 2 ^{h - h) = 

I2UJ2 + ^3^i(h — h) — (3.33) 
I 3 u 3 + LO X Lo 2 {h ~ h) = 

These are Euler's Equations. 

We can extend this analysis to include a torque r. The equation of motion becomes 
L = r and we can again expand in the body frame along the principal axes to derive 
Euler's equations (3.33), now with the components of the torque on the RHS. 

3.4 Free Tops 

"To those who study the progress of exact science, the common spinning-top 
is a symbol of the labours and the perplexities of men." 

James Clerk Maxwell, no less 

In this section, we'll analyse the motion of free rotating bod- 
ies (known as free tops) using Euler's equation. 

We start with a trivial example: the sphere. For this object, 
h — h — h which means that the angular velocity u) is parallel 
to the angular momentum L. Indeed, Euler's equations tell us 
that uj a is a constant in this case and the sphere continues to spin 
around the same axis you start it on. To find a more interesting 
case, we need to look at the next simplest object. 

3.4.1 The Symmetric Top Figure 29: 

The symmetric top is an object with I x = I 2 ^ I 3 . A typical example is drawn in figure 
29. Euler's equations become 

hu-L = U 2 U2,{h ~ h) 

I2UJ2 = -wiu 3 (Ii - I 3 ) (3.34) 
I 3 u 3 = 
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Figure 30: The precession of the spin: the direction of precession depends on whether the 
object is short and fat (I3 > h) or tall and skinny (I 3 < I\) 



So, in this case, we see that uj 3 , which is the spin about the symmetric axis, is a 
constant of motion. In contrast, the spins about the other two axes are time dependant 
and satisfy 



U02 = — 



(3.35) 



where 



n = cu 3 (h - h)/h 



(3.36) 



is a constant. These equations are solved by 



(u>i,u>2) — ou (smQt,cosQt) 



(3.37) 



for any constant ujq. This means that, in the body frame, the 
direction of the spin is not constant: it precesses about the e 3 axis 
with frequency Q. The direction of the spin depends on the sign on 
Q or, in other words, whether I\ > I 3 or I\ < I 3 . This is drawn in 
figure 30. 

In an inertial frame, this precession of the spin looks like a wobble. 
To see this, recall that L has a fixed direction. Since both oo 3 and 
L3 are constant in time, the e 3 axis must stay at a fixed angle with 
respect to the L and u). It rotates about the L axis as shown in figure 
31. We'll examine this wobble more in the next section. 




Figure 31: 
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3.4.2 Example: The Earth's Wobble 

The spin of the earth causes it to bulge at the equator so it is no longer a sphere but can 
be treated as a symmetric top. It is an oblate ellipsoid, with 7 3 > I ± , and is spherical 
to roughly 1 part in 300, meaning 

^— ^ » -— (3.38) 
h 300 1 ; 

Of course, we know the magnitude of the spin cj 3 : it is uo% = (1 day) -1 . This information 

is enough to calculate the frequency of the earth's wobble; from (3.36), it should be 

1 



earth 



300 



day 



-i 



(3.39) 




Equator 



This calculation was first performed by Euler in 1749 who pre- 
dicted that the Earth completes a wobble every 300 days. De- 
spite many searches, this effect wasn't detected until 1891 when 
Chandler re-analysed the data and saw a wobble with a period 
of 427 days. It is now known as the Chandler wobble. It is 
very small! The angular velocity u> intercepts the surface of 
the earth approximately 10 metres from the North pole and 
precesses around it. More recent measurements place the fre- 
quency at 435 days, with the discrepancy between the predicted 
300 days and observed 435 days due to the fact that the Earth 
is not a rigid body, but is flexible because of tidal effects. Less 
well understood is why these same tidal effects haven't caused 

the wobble to dampen and disappear completely. There are various theories about 
what keeps the wobble alive, from earthquakes to fluctuating pressure at the bottom 
of the ocean. 

3.4.3 The Asymmetric Top: Stability 

The most general body has no symmetries and h ^ h ^ h ^ h- The rotational 
motion is more complicated but there is a simple result that we will describe here. 
Consider the case where the spin is completely about one of the principal axes, say e^ 
i.e. 



Figure 32: 



u 2 = uj 3 = 



(3.40) 



This solves Euler's equations (3.33). The question we want to ask is: what happens if 
the spin varies slightly from this direction? To answer this, consider small perturbations 
about the spin 



^1 = ^ + 77! 



uj 2 = r)2 



^3 = ??3 



(3.41) 
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where r] a , a — 1, 2, 3 are all taken to be small. Substituting this into Euler's equations 
and ignoring terms of order rf and higher, we have 



hm = o 

I 2 fj 2 = ^ 3 (/ 3 - h) 

I 3 rj 3 = Qrj 2 (Ii - h) 



(3.42) 
(3.43) 



We substitute the third equation into the second to find an equation for just one of the 
perturbations, say r) 2 , 



The fate of the small perturbation depends on the sign of the quantity A. We have two 
possibilities 

• A < 0: In this case, the disturbance will oscillate around the constant motion. 

• A > 0: In this case, the disturbance will grow exponentially. 
Examining the definition of A, we find that the motion is unstable if 



with all other motions stable. In other words, a body will rotate stably about the axis 
with the largest or the smallest moment of inertia, but not about the intermediate axis. 
Pick up a tennis racket and try it for yourself! 

3.4.4 The Asymmetric Top: Poinsot Construction 

The analytic solution for the general motion of an asymmetric top is rather complicated, 
involving Jacobian elliptic functions. But there's a nice geometrical way of viewing the 
motion due to Poinsot. 

We start by working in the body frame. There are two constants of motion: the 
kinetic energy T and the magnitude of the angular momentum L 2 . In terms of the 
angular velocity, they are 



hm = -r(h ~ h){h ~ h)m = Ar] 2 



(3.44) 



I 2 < h< h or I 3 < h < I 2 



(3.45) 



2T = hul + I 2 uj\ + hul 



(3.46) 
(3.47) 
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Each of these equations defines an ellipsoid in uj space. The motion of the vector uj is 
constrained to lie on the intersection of these two ellipsoids. The first of these ellipsoids, 
defined by 



is known as the inertia ellipsoid (or, sometimes, the inertia quadric). If we fix the 
kinetic energy, we can think of this abstract ellipsoid as embedded within the object, 
rotating with it. 



circles around the u)\ and 003 axes, but two 

lines passing through the o>2 axis. For fixed T and L 2 , the vector uj moves along one 
of the intersection lines. This provides a pictorial demonstration of the fact we learnt 
in the previous subsection: an object will spin in a stable manner around the principal 
axes with the smallest and largest moments of inertia, but not around the intermediate 
axis. The path that uj traces on the inertia ellipsoid is known as the polhode curve. We 
see from the figure that the polhode curves are always closed, and motion in the body 
frame is periodic. 

So much for the body frame. What does all this look like in the space frame? The 
vector L is a constant of motion. Since the kinetic energy 2T = L • uj is also constant, 
we learn that uj must lie in a fixed plane perpendicular to L. This is known as the 
invariable plane. The inertia ellipsoid touches the invariable plane at the point defined 
by the angular velocity vector uj. Moreover, the invariable plane is always tangent to 
the inertial ellipsoid at the point uj. To see this, note that the angular momentum can 
be written as 




(3.48) 



The inertia ellipsoid is drawn in figure 
33, where we've chosen 1\ > I 2 > I3 so that 



A 



the major axis is 003 and the minor axis is uj\. 
The lines drawn on the figure are the inter- 
section of the inertia ellipsoid with the other 
ellipsoid, defined by (3.47), for various val- 
ues of L 2 . Since this has the same major and 
minor axes as the inertia ellipsoid (because 
If > If > If), the intersection lines are small 




Figure 33: 



(3.49) 



where the gradient operator is in uj space, i.e. V w = (d/dui,d/duj2,d/du)3). But 
recall that the inertia ellipsoid is defined as a level surface of T, so equation (3.49) tells 




Figure 34: The inertia ellipsoid rolling around on the invariable plane, with the polhode 
and herpolhode curves drawn for a fixed time period. 

us that the angular momentum L is always perpendicular to the ellipsoid. This, in 
turn, ensures that the invariable plane is always tangent to the ellipsoid. In summary, 
the angular velocity traces out two curves: one on the inertia ellipsoid, known as the 
polhode curve, and another on the invariable plane, known as the herpolhode curve. 
The body moves as if it is embedded within the inertia ellipsoid, which rolls around the 
invariable plane without slipping, with the center of the ellipsoid a constant distance 
from the plane. The motion is shown in figure 34. Unlike the polhode curve, the 
herpolhode curve does not necessarily close. 

An Example: The Asteroid Toutatis 

Astronomical objects are usually symmetric, but there's an important exception wan- 
dering around our solar system, depicted in figure 2 35. This is the asteroid Toutatis. In 
September 2004 it passed the earth at a distance of about four times that to the moon. 
This is (hopefully!) the closest any asteroid will come for the next 60 years. The orbit 
of Toutatis is thought to be chaotic, which could potentially be bad news for Earth 
a few centuries from now. As you can see from the picture, its tumbling motion is 
complicated. It is aperiodic. The pictures show the asteroid at intervals of a day. The 
angular momentum vector L remains fixed and vertical throughout the motion. The 
angular velocity uj traces out the herpolhode curve in the horizontal plane, perpendic- 
ular to L. The angular momentum vector u> also traces out a curve over the asteroid's 

2 This picture was created by Scott Hudson of Washington State University and was taken from 
http://www.solarviews.com/eng/toutatis.htm where you can find many interesting facts about the 
asteroid. 



- 60 - 



Figure 35: By Toutatis! The three principal axes are shown in red, green and blue (without 
arrows). The angular momentum L is the vertical, purple arrow. The angular velocity u is 
the circled, yellow arrow. 

surface: this is the polhode curve. It has a period of 5.4 days which you can observe 
by noting that u) has roughly the same orientation relative to the principal axes every 
five to six days. 

3.5 Euler's Angles 

So far we've managed to make quite a lot of progress working just with the angular 
velocity uo a and we haven't needed to introduce an explicit parametrization of the 
configuration space C . But to make further progress we're going to need to do this. 
We will use a choice due to Euler which often leads to simple solutions. 
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Figure 36: The rotation from space frame {e a } to body frame {e a }. 

A general rotation of a set of axis is shown in Figure 36. We'd like to construct a 
way of parameterizing such a rotation. The way to do this was first described by Euler: 

Euler's Theorem: 

An arbitrary rotation may be expressed as the product of 3 successive rotations about 
3 (in general) different axes. 

Proof: Let {e a } be space frame axes. Let {e a } be body frame axes. We want to find 
the rotation R so that e a = R a b^b- We can accomplish this in three steps 



l e 4 — ► {ej — ► {ej — ► {e a \ 



(3.50) 



Let's look at these step in turn. 




Figure 37: Step 1: Rotate around the space-frame axis §3. 

Step 1: Rotate by <fi about the e 3 axis. So e' a = R^,{(j))ab^b with 

^ cos0 sin0 o\ 

RM = 



— sin <f> cos 4> 
1 



(3.51) 
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Figure 38: Step 2: Rotate around the new axis axis e[. 



Step 2: Rotate by 9 about the new axis e[. This axis e' : is sometimes called the 
"line of nodes". We write e" = Ri(9)e' b with 

/l \ 



Ri{9) 



cos 9 sin 9 
— sin 9 cos 9 



This is shown in Figure 38 





A 



(3.52) 



Figure 39: Step 3: Rotate around the latest axis e(j. 



Step 3: Rotate by ip about the new new axis e^' so e a = R 3 (ip) ab e' b ' with 

^ cos^ sinip o\ 



sin if; cos tp 
1 



This is shown in Figure 39. 



(3.53) 



- 63 - 



So putting it all together, we have 

iU&M) = [RsWM8)Rs(<f>)L b (3-54) 

□ 

The angles 0, 9, ip are the Euler angles. If we write out the matrix i?(0, 9, ip) longhand, 
it reads 

^ cos ip cos — cos 9 sin sin ^ sin cos ^ + cos 9 sin ^ cos sin # sin ^ \ 
R — — cos sin *0 — cos cos ^ sin — sin -0 sin + cos cos -0 cos sin cos -0 
^ sin sin — sin 9 cos cos 9 j 

Note: Recall that we may expand a vector r either in the body frame r = r a e a , or in 
the space frame r = r a e a . The above rotations can be equally well expressed in terms 
of the coordinates r a rather than the basis {e a }: we have f b = r a R ab . Be aware that 
some books choose to describe the Euler angles in terms of the coordinates r a which 
they write in vector form. In some conventions this can lead to an apparent reversal in 
the ordering of the three rotation matrices. 

3.5.1 Leonhard Euler (1707-1783) 

As is clear from the section headings, the main man for this chapter is Euler, by far the 
most prolific mathematician of all time. As well as developing the dynamics of rotations, 
he made huge contributions to the fields of number theory, geometry, topology, analysis 
and fluid dynamics. For example, the lovely equation e t6 = cos 9 + isin 9 is due to 
Euler. In 1744 he was the first to correctly present a limited example of the calculus of 
variations (which we saw in section 2.1) although he generously gives credit to a rather 
botched attempt by his friend Maupertuis in the same year. Euler also invented much 
of the modern notation of mathematics: f(x) for a function; e for exponential; it for, 
well, 7i and so on. 

Euler was born in Basel, Switzerland and held positions in St Petersburg, Berlin and, 
after falling out with Frederick the Great, St Petersburg again. He was pretty much 
absorbed with mathematics day and night. Upon losing the sight in his right eye in 
his twenties he responded with: "Now I will have less distraction" . Even when he went 
completely blind later in life, it didn't slow him down much as he went on to produce 
over half of his total work. The St Petersburg Academy of Science continued to publish 
his work for a full 50 years after his death. 

3.5.2 Angular Velocity 

There's a simple expression for the instantaneous angular velocity u) in terms of Euler 
angles. To derive this, we could simply plug (3.54) into the definition of angular velocity 
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(3.5). But this is tedious, and a little bit of thought about what this means physically 
will get us there quicker. Consider the motion of a rigid body in an infinitesimal time 
dt during which 

(ip,9,<f>) -> (V> + #,0 + dM + #) (3-55) 
From the definition of the Euler angles, the angular velocity must be of the form 

w = <j)e 3 + 9e[ +^e 3 (3.56) 
But we can express the first two vectors in terms of the body frame. They are 

e 3 = sin 9 sin ipei + sin 9 cos vpe 2 + cos 9 e 3 

e[ = cos ip e i — sin ^ e 2 (3.57) 
from which we can express w in terms of the Euler angles in the body frame axis 

u) = [(p sin 9 sin tp + 9 cos ip]ei + [<p sin 9 cos tp — 9 sin T/>]e 2 + + cos #]e 3 (3.58) 
By playing a similar game, we can also express ui in the space frame axis. 
3.5.3 The Free Symmetric Top Revisited 

In section 3.4 we studied the free symmetric top working in the body frame and found a 
constant spin u; 3 while, as shown in equation (3.37), lo\ and 0J2 precess with frequency 

n = uj 3 ^ h 7 /3) (3.59) 

h 

But what does this look like in the space framel Now that we have parametrised 
motion in the space frame in terms of Euler angles, we can answer this question. This 
are simplest if we choose the angular momentum L to lie along the e 3 space-axis. Then, 
since we have already seen that e 3 sits at a fixed angle to L, from the figure we see 
that 9 = 0. Now we could either use the equations (3.58) or, alternatively, just look at 
figure 40, to see that we should identify Q — ip. 

But we know from (3.58) that the expression for u 3 (which, remember, is the compo- 
nent of u> in the body frame) in terms of Euler angles is u; 3 = ip + <j) cos 9 so, substituting 
for Q = ip, we find the precession frequency 



Figure 40: Euler angles for the free symmetric top when L coincides with £3 

An Example: The Wobbling Plate 

The physicist Richard Feynman tells the following story: 

"J was in the cafeteria and some guy, fooling around, throws a plate in the 
air. As the plate went up in the air I saw it wobble, and I noticed the red 
medallion of Cornell on the plate going around. It was pretty obvious to me 
that the medallion went around faster than the wobbling. 

I had nothing to do, so I start figuring out the motion of the rotating plate. 
I discover that when the angle is very slight, the medallion rotates twice as 
fast as the wobble rate - two to one. It came out of a complicated equation! 

I went on to work out equations for wobbles. Then I thought about how the 
electron orbits start to move in relativity. Then there's the Dirac equation 
in electrodynamics. And then quantum electrodynamics. And before I knew 
it.... the whole business that I got the Nobel prize for came from that piddling 
around with the wobbling plate." 

Feynman was right about quantum electrodynamics. But what about the plate? 
We can look at this easily using what we've learnt. The spin of the plate is u 3 , while 
the precession, or wobble, rate is which is given in (3.60). To calculate this, we need 
the moments of inertia for a plate. But we figured this out for the disc in Section 3.2 
where we found that I3 = 2I±. We can use this to see that tp = —U3 for this example 
and so, for slight angles 9, have 

0^-2^ (3.61) 
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Or, in other words, the wobble rate is twice as fast as the spin 
of the plate. It's the opposite to how Feynman remembers! 

There is another elegant and simple method you can use to see 
that Feynman was wrong: you can pick up a plate and throw it. It's 
hard to see that the wobble to spin ratio is exactly two. But it's easy 
to see that it wobbles faster than it spins. 




Figure 41: 

3.6 The Heavy Symmetric Top 




Figure 42: The heavy top with its Euler angles 

The "heavy" in the title of this section means that the top is acted upon by gravity. 
We'll deal only with a symmetric top, pinned at a point P which is a distance I from 
the centre of mass. This system is drawn in the figure. The principal axes are ei, e2 
and e 3 and we have I\ — I 2 - From what we have learnt so far, it is easy to write down 
the Lagrangian: 

L = + ) + \hw\ - Mgl cos 9 

= \h(6 2 + sin 2 6<j) 2 ) + §/ 3 (V> + cos 6 0) 2 - Mgl cos 6 (3.62) 

A quick examination of this equation tells us that both ip and are ignorable coordi- 
nates. This gives us the constants of motion and p^, where 

p^ = I 3 (ip + cos 9 4>)= I 3 u 3 (3.63) 
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This is the angular momentum about the symmetry axis e3 of the top. The angular 
velocity 0J3 about this axis is simply called the spin of the top and, as for the free 
symmetric top, it is a constant. The other constant of motion is 

p<p = h sin 2 9 + I 3 cos 9 (ip + cos 9) (3.64) 

As well as these two conjugate momenta, the total energy E is also conserved 

E = T + V = \h (9 2 + 2 sin 2 9) + |/ 3 cj 2 + Mgl cos 9 (3.65) 

To simplify these equations, let's define the two constants 



Then we can write 



and 



and b = ^ (3.66) 



X b-acos9 

= -T7T- ( 3 - 67 ) 

sin 



ha (b - a cos 9) cos 9 

^ = -j — a (3-68) 

h sin 9 



So if we can solve 9 = 9{t) somehow, then we can always integrate these two equations 
to get <f)(t) and ip(t). But first we have to figure out what 9 is doing. To do this, let's 
define the "reduced energy" E' = E — ^hwl- Then, since E and u 3 are constant, so is 
E'. We have 

E' = l -h9 2 + V eS (9) (3.69) 

where the effective potential is 

VM<»= h(b : a ?: er ' +M alm e (3.70) 
2 sin 9 

So we've succeeded in getting an equation (3.69) purely in terms of 9. To simplify the 
analysis, let's define the new coordinate 

w = cos# (3.71) 

Clearly —1 < u < 1. We'll also define two further constants to help put the equations 
in the most concise form 

a = — and (3 = (3.72) 
h h 
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With all these redefinitions, the equations of motion (3.67), (3.68) and (3.69) can be 
written as 



u 2 = (1 - u 2 )(a - (3u) - (b - au) 2 = f(u) 

■ b — au 

<f) = 



1-u 2 
Iia u(b — au) 



1-u 2 



(3.73) 
(3.74) 

(3.75) 



We could take the square root of equation (3.73) and integrate to reduce the problem 
to quadrature. The resulting integral is known as an "elliptic integral". But, rather 
than doing this, there's a better way to understand the physics qualitatively. 

Note that the function f(u) defined in (3.73) is a cubic polynomial that behaves as 



/(«) 



+oo as u — > oo 



-oo as u 



(3.76) 



-oo 



and /(±1) = —(b =F a) 2 < 0. So if we plot the function f(u), it looks like figure 43 




Figure 43: 



The physical range is u 2 = f(u) > and — 1 < u < 1 so we find that, like in the 
spherical pendulum and central force problem, the system is confined to lie between 
the two roots of f(u). 

There are three possibilities for the motion depending on the sign of <j> at the two 
roots u — ui and u = u<i as determined by (3.74). These are 

• > at both u — U\ and u = u 2 
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Motion 1) 



Motion 2) 



Motion 3) 




Figure 44: The three different types of motion depend on the direction of precession at the 
extremal points. 

• (f> > at u — ui, but < at u = u 2 

• <fi > at u = ui and = at u = u 2 

The different paths of the top corresponding to these three possibilities are shown in 
figure 44. Motion in <p is called precession while motion in 9 is known as nutation. 

3.6.1 Letting the Top go 

The last of these three motions is not as unlikely as it may first appear. Suppose 
we spin the top and let it go at some angle 9. What happens? We have the initial 
conditions 

t=o = O /(«t=o) = 

=>• u t=0 = u 2 
and <p t =o = =3- b — au t=0 = 

«t=o = - (3.77) 
a 

Remember also that the quantity 

Ptl, = h<j) sin 2 9 + I 3 u; 3 cos 9 = I 3 lu 3 cos 9 t=0 (3.78) 

is a constant of motion. We now have enough information to figure out the qualitative 
motion of the top. Firstly, it starts to fall under the influence of gravity, so 9 increases. 
But as the top falls, <fi must turn and increase in order to keep p<f, constant. Moreover, 
we also see that the direction of the precession <f> must be in the same direction as the 
spin u 3 itself. What we get is motion of the third kind. 
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3.6.2 Uniform Precession 



(3.79) 



Can we make the top precess with bobbing up and down? i.e. with 9 = and (ft 
constant. For this to happen, we would need the function f(u) to have a single root uq 
lying in the physical range — 1 < uq < +1. This root must satisfy 

f(u ) = (1 - ul)(a - (3u Q ) - (b - au ) 2 = 
and f'(uo) = = —2u (a — (5uq) — (3(1 — ufy + 2a{b — au ) = 

Combining these, we find \fi = acp — <p 2 uo. Substitut- m 
ing the definitions I±a = I3U3 and (5 = 1MgljI\ into this 
expression, we find 

Mgl = 0(/ 3 ^3 - h4> cos Q ) (3.80) 

The interpretation of this equation is as follows: for a fixed 
value of U3 (the spin of the top) and 9 (the angle at which 
you let it go), we need to give exactly the right push <fi to Figure 45: 

make the top spin without bobbing. In fact, since equation (3.80) is quadratic in 0, 
there are two frequencies with which the top can precess without bobbing. 

Of course, these "slow" and "fast" precessions only 
exist if equation (3.80) has any solutions at all. Since it is 
quadratic, this is not guaranteed, but requires 




u 3 > — \/ Mglli cos 9 
I?. 



(3.81) 



So we see that, for a given #o> the top has to be spinning fast 
enough in order to have uniform solutions. What happens 
if it's not spinning fast enough? Well, the top falls over! 

3.6.3 The Sleeping Top 

Suppose we start the top spinning in an upright position, with 

9 = 9 = 




Figure 46: 



(3.82) 



When it spins upright, it is called a sleeping top. The question we want to answer is: 
will it stay there? Or will it fall over? From (3.73), we see that the function f(u) must 
have a root at 9 = 0, or u = +1: f(l) = 0, From the definitions (3.66) and (3.72), we 
can check that a = b and a = j3 in this situation and f(u) actually has a double pole 
at u = +1, 

f(u) = (1 -u) 2 (a(l + u) -a 2 ) (3.83) 
The second root of f(u) is at u 2 = a 2 /a — 1. There are two possibilities 
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Figure 47: The function f(u) for the stable and unstable sleeping top. 

1: u 2 > 1 or cj| > AIiMgl/I^. In this case, the graph of f(u) is drawn in first in 
figure 47. This motion is stable: if we perturb the initial conditions slightly, we will 
perturb the function f(u) slightly, but the physical condition that we must restrict to 
the regime f{u) > means that the motion will continue to be trapped near u = 1 

2: u 2 < 1 or cjf < 4i 1 Mg///f . In this case, the function f(u) looks like the second 
figure of 47. Now the top is unstable; slight changes in the initial condition allow a 
large excursion. 

In practice, the top spins upright until it is slowed by friction to I 3 u 3 = 2^/IiMgl, 
at which point it starts to fall and precess. 

3.6.4 The Precession of the Equinox 

The Euler angles for the earth are drawn in figure 48. The earth spins at an angle of 
9 = 23.5° to the plane of its orbit around the sun (known as the plane of the elliptic). 
The spin of the earth is ip — (day) -1 . This causes the earth to bulge at the equator 
so it is no longer a sphere, but rather a symmetric top. In turn, this allows the moon 
and sun to exert a torque on the earth which produces a precession <fi. Physically this 
means that the direction in which the north pole points traces a circle in the sky and 
what we currently call the "pole star" will no longer be in several thousand years time. 
It turns out that this precession is "retrograde" i.e. opposite to the direction of the 
spin. 

One can calculate the precession <fi of the earth due to the moon and sun using the 
techniques described in the chapter. But the calculation is rather long and we won't go 
over it in this course (see the book by Hand and Finch if you're interested). Instead, 
we will use a different technique to calculate the precession of the earth: astrology! 3 

3 I learnt about this fact from John Baez' website where you can find lots of well written explanations 
of curiosities in mathematical physics: http://math.ucr.edu/home/baez/README.html. 
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June 23 




Eailh 




Sun 



Dec 22 




Sep 23 



Figure 48: The precession of the earth. 



To compute the precession of the earth, the first fact we need to know is that Jesus 
was born in the age of Pisces. This doesn't mean that Jesus looked up Pisces in his 
daily horoscope (while scholars are divided over the exact date of his birth, he seems to 
exhibit many traits of a typical Capricorn) but rather refers to the patch of the sky in 
which the sun appears during the first day of spring. Known in astronomical terms as 
the "vernal equinox", this day of the year is defined by the property that the sun sits 
directly above the equator at midday. As the earth precesses, this event takes place 
at a slightly different point in its orbit each year, with a slightly different backdrop of 
st result. The astrological age is defined to be the background constellation in 

which the sun rises during vernal equinox. 

It is easy to remember that Jesus was born in the age of Pisces since the fish was 
used as an early symbol for Christianity. The next fact that we need to know is 
that we're currently entering the age of Aquarius (which anyone who has seen the 
musical Hair will know). So we managed to travel backwards one house of the zodiac 
in 2,000 years. We've got to make it around 12 in total, giving us a precession time of 
2, 000 x 12 = 24, 000 years. The actual value of the precession is 25, 700 years. Our 
calculation is pretty close considering the method! 

The earth also undergoes other motion. The value of 9 varies from 22.1° to 24.5° 
over a period of 41, 000 years, mostly due to the effects of the other planets. These also 
affect the eccentricity of the orbit over a period of 105,000 years. 

3.7 The Motion of Deformable Bodies 

Take a lively cat. (Not one that's half dead like Schrodinger's). Hold it upside down and 
drop it. The cat will twist its body and land sprightly on its feet. Yet it doesn't do this 
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by pushing against anything and its angular momentum is zero throughout. If the cat 
were rigid, such motion would be impossible since a change in orientation for a rigid 
body necessarily requires non-vanishing angular momentum. But the cat isn't rigid 
(indeed, it can be checked that dead cats are unable to perform this feat) and bodies 
that can deform are able to reorient themselves without violating the conservation of 
angular momentum. In this section we'll describe some of the beautiful mathematics 
that lies behind this. I should warn you that this material is somewhat more advanced 
than the motion of rigid bodies. The theory described below was first developed in the 
late 1980s in order to understand how micro-organisms swim 4 . 

3.7.1 Kinematics 

We first need to describe the configuration space C of a deformable body. We factor 
out translations by insisting that all bodies have the same center of mass. Then the 
configuration space C is the space of all shapes with some orientation. 

Rotations act naturally on the space C (they simply rotate each shape). This allows 
us to define the smaller shape space C so that any two configurations in C which are 
related by a rotation are identified in C. In other words, any two objects that have the 
same shape, but different orientation, are described by different points in C, but the 
same point in C. Mathematically, we say C = C/SO(3). 

We can describe this in more detail for a body consisting of N point masses, each with 
position rj. Unlike in section 3.1, we do not require that the distances between particles 
are fixed, i.e. |rj — Tj\ 7^ constant. (However, there may still be some restrictions on 
the i^). The configuration space C is the space of all possible configurations rj. For 
each different shape in C, we pick a representative f j with some, fixed orientation. It 
doesn't matter what representative we choose — just as long as we pick one. These 
variables fj are coordinates on the space shape C. For each r^ e C, we can always find 
a rotation matrix R e 5*0(3) such that 

v i = Ri l (3.84) 

As in section 3.1, we can always this to continuous bodies. In this case, the configuration 
space C and the shape space C may be infinite dimensional. Examples of different shapes 
for a continuously deformable body are shown in figure 49. 



4 See A. Shapcrc and F. Wilczek, " Geometry of Self- Propulsion at Low Reynolds Number" , J. Fluid 
Mech. 198 557 (1989) . This is the same Frank Wilczek who won the 2004 Nobel prize for his work 
on quark interactions. 
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Figure 49: Three possible shapes of a deformable object. 



We want to understand how an object rotates as it changes shape keeping its 
angular momentum fixed (for example, keeping L = throughout). The first thing to 
note is that we can't really talk about the rotation between objects of different shapes. 
(How would you say that the the third object in figure 49 is rotated with respect to the 
first or the second?). Instead, we should think of an object moving through a sequence 
of shapes before returning to its initial shape. We can then ask if there's been a net 
rotation. As the object moves through its sequence of shapes, the motion is described 
by a time dependent fj(t), while the corresponding change through the configuration 
space is 

Ti(t) = R(t) r(t) (3.85) 

where the 3x3 rotation matrix R(t) describes the necessary rotation to go from our 
fixed orientation of the shape f to the true orientation. As in section 3.1.1, we can define 
the 3x3 anti-symmetric matrix that describes the instantaneous angular velocity of 
the object. In fact, it will for once prove more useful to work with the "connective 
angular velocity" defined around equation (3.10) 

« = (3.86) 

This angular velocity is non-zero due to the changing shape of the object, rather than 
the rigid rotation that we saw before. Let's do a quick change of notation and write 
coordinates on the shape space C as x A , with A = 1,...,3N instead of in vector 
notation f with i — 1, . . . , N. Then, since Q is linear in time derivatives, we can write 

n = n A (x)x A (3.87) 

The component Qa(x) is the 3x3 angular velocity matrix induced if the shape changes 
from x A to x A + 5x A . It is independent of time: all the time dependence sits in 
the x A factor which tells us how the shape is changing. The upshot is that for each 
shape x G C, we have a 3 x 3 anti-symmetric matrix Q A associated to each of the 
A = 1, . . . , 3N directions in which the shape can change. 
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However, there is an ambiguity in defining the angular velocity f2. This comes about 
because of our arbitrary choice of reference orientation when we picked a representative 
f j e C for each shape. We could quite easily have picked a different orientation, 

f . _> S(x A ) ii (3.88) 

where S(x a ) is a rotation that, as the notation suggests, can vary for each shape x A . If 
we pick this new set of representative orientations, then the rotation matrix R defined 
in (3.85) changes: R(t) — > R(t) S~ 1 (x A ). Equation (3.86) then tells us that the angular 
velocity also change as 

Qa^SQaS-' + S^ (3.89) 

This ambiguity is related to the fact that we can't define the meaning of rotation 
between two different shapes. Nonetheless, we will see shortly that when we come to 
compute the net rotation of the same shape, this ambiguity will disappear, as it must. 
Objects such as Qa which suffer an ambiguity of form (3.89) are extremely important 
in modern physics and geometry. They are known as non-abelian gauge potentials to 
physicists, or as connections to mathematicians. 

3.7.2 Dynamics 

So far we've learnt how to describe the angular velocity Q of a deformable object. The 
next step is to see how to calculate Q. We'll now show that, up to the ambiguity 
described in (3.89), the angular velocity Q is specified by the requirement that the 
angular momentum L of the object is zero. 

L = ^m ! r i xr i 

i 

= ^m f [(iffi) x (Rii) + (Rrt) x (i?f,)] = (3.90) 

i 

In components this reads 

L a = eabc^mi ^R bd R ce (r i ) d (r i ) e + i2 M ^ ce (f i ) d (f i ) c ] = (3.91) 

i 

The vanishing L = is enough information to determine the following result: 
Claim: The 3x3 angular velocity matrix Vl ab = R~ l R cb is given by 

&<ab = tabc I c d Ld (3.92) 
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where / is the instantaneous inertia tensor of the shape described by fj, 

Lb = ^^((ii ■ Yi)5 ab - (f i ) (f i ) 6 ) (3.93) 

i 

and L a is the apparent angular momentum 

~L a = e a bc ^2 mi ( f <)ft( f <)c ( 3 - 94 ) 



Proof: We start by multiplying L a by e a f g . We need to use the fact that if we 
multiply two e-symbols, we have e abc e a fg = (Sbf^cg ~ 8bg8 c f)- Then 

e afgL a = ^2 m i [Rfd,Rge{ii)d{ii)e ~ RgdRfe{ii)d{ii)e 

i 

-RgdRfe(Ti)d(Ti)e + RfdRge(ii)d(ii)e = (3.95) 

Now multiply by Rf b R gc . Since R is orthogonal, we known that Rf b R/d = hd which, 
after contracting a bunch of indices, gives us 

RfbRg C e a fgL a = "^2 m i [( r i)&( r i)c - (fi)c(fi) 6 - n M (f i ) c (f i ) d + n cd (f i ) 6 (f i ) d ] = 

This is almost in the form that we want, but the indices aren't quite contracted in the 
right manner to reproduce (3.92). One can try to play around to get the indices working 
right, but at this stage it's just as easy to expand out the components explicitly. For 
example, we can look at 

L l = ^m i [(f i ) 2 (f i ) 3 - ((f i ) 3 (f i ) 2 ] 

i 

= ^ m i [ fi 2l(fi) 3 (fj)l + ^23(^)3(^)3 - fi 3 l(fj) 2 (fi)l - ^32(^)2(^)2] 

i 

= ^11^23 + ^12^31 + ^13^12 = \ tabchaQbc (3.96) 

where the first equality is the definition of L ± , while the second equality uses our result 
above, and the third equality uses the definition of / given in (3.93). There are two 
similar equations, which are summarised in the formula 

La — \ tbcdlab^cd (3.97) 

Multiplying both sides by J -1 gives us precisely the claimed result (3.92). This con- 
cludes the proof. □. 
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To summarise: a system with no angular momentum that can twist and turn and 
change its shape has an angular velocity (3.92) where f i(t) is the path it chooses to take 
through the space of shapes. This is a nice formula. But what do we do with it? We 
want to compute the net rotation R as the body moves through a sequence of shapes 
and returns to its starting point at a time T later. This is given by solving (3.86) for 
R. The way to do this was described in section 3.1.2. We use path ordered exponentials 



The path ordering symbol P puts all matrices evaluated at later times to the right. 
(This differs from the ordering in section 3.1.2 where we put later matrices to the left. 



exponent as an integral around a closed path in shape space. Here time has dropped 
out. This tells us an important fact: it doesn't matter how quickly we perform the 
change of shapes — the net rotation of the object will be the same. 

In particle physics language, the integral in (3.98) is called a "Wilson loop". We can 
see how the rotation fares under the ambiguity (3.87). After some algebra, you can 
find that the net rotation R of an object with shape x A is changed by 



This is as it should be: the S' 1 takes the shape back to our initial choice of standard 
orientation; the matrix R is the rotation due to the change in shape; finally S puts us 
back to the new, standard orientation. So we see that even though the definition of the 
angular velocity is plagued with ambiguity, when we come to ask physically meaningful 
questions - - such as how much has a shape rotated — the ambiguity disappears. 
However, if we ask nonsensical questions — such as the rotation between two different 
shapes — then the ambiguity looms large. In this manner, the theory contains a rather 
astonishing new ingredient: it lets us know what are the sensible questions to ask! 
Quantities for which the ambiguity (3.87) vanishes are called gauge invariant. 

In general, it's quite hard to explicitly compute the integral (3.98). One case where 
it is possible is for infinitesimal changes of shape. Suppose we start with a particular 
shape x A , and move infmitesimally in a loop in shape space: 




(3.98) 




R^ S(x A )RS(x A )- 



i 



(3.99) 



x A (t) =x° A + a A (t) 



(3.100) 



Then we can Taylor expand our angular velocity components, 



n A (x(t)) = n A (x°) + 



dtt A 



dx B 



a B (3.101) 



,.0 



Expanding out the rotation matrix (3.98) and taking care with the ordering, one can 
show that 

R=l + \ f ab j a A a B dt + 0(a 3 ) 

= 1 + l - J F AB dA AB + 0(a 3 ) (3.102) 

where F AB is anti-symmetric in the shape space indices A and B, and is a 3 x 3 matrix 
(the a,b — 1,2,3 indices have been suppressed) given by 

^ = 0-^ + ^] < 3 - 103 > 

It is known as the field strength to physicists (or the curvature to mathematicians). It is 
evaluated on the initial shape x A . The second equality in (3.102) gives the infinitesimal 
rotation as the integral of the field strength over the area traversed in shape space. This 
field strength contains all the information one needs to know about the infinitesimal 
rotations of objects induced by changing their shape. 

One of the nicest things about the formalism described above is that it mirrors very 
closely the mathematics needed to describe the fundamental laws of nature, such as 
the strong and weak nuclear forces and gravity. They are all described by "non-abelian 
gauge theories", with an object known as the gauge potential (analogous to Q A ) and 
an associated field strength. 
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